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Abstract. The purpose of this paper is describe Lagrangian Mechanics for 
constrained systems on Lie algebroids, a natural framework which covers a 
wide range of situations (systems on Lie groups, quotients by the action of 
a Lie group, standard tangent bundles...). In particular, we are interested in 
two cases: singular Lagrangian systems and vakonomic mechanics (varia- 
tional constrained mechanics) . Several examples illustrate the interest of these 
developments. 
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1. Introduction 

There is a vast literature around the Lagrangian formalism in mechanics justified 
by the central role played by these systems in the foundations of modern mathema- 
tics and physics. In some interesting systems some problems often arise due to their 
singular nature that give rise to the presence of constraints manifesting that the 
evolution problem is not well posed {internal constraints). Constraints can also 
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manifest a priori restrictions on the states of the system which are often imposed 
either by physical arguments or by external conditions {external constraints). 
Both cases are of considerable importance. 

Systems with internal constraints are quite interesting since many dynamical 
systems are given in terms of non-symplectic forms instead of the more habitual 
symplectic ones. The more frequent case appears in the Lagrangian formalism of 
singular mechanical systems which are a commonplace in many physical theories 
(as in Yang-Mills theories, gravitation, etc). Also, in systems that appears as a limit 
(as in Chern-Simons lagrangians); for instance consider the following lagrangian 

L = ^m,iq'f + eMq)q'-V{q) 

in the limit nii for some i. In other cases, it is necessary work with a new 
singular lagrangian in a extended space since the original lagrangian is "ill defined" 
(only locally defined on the original space) , as it happens for the electron monopole 
system (see [16]). 

Another motivation for the present work is the study of lagrangian systems sub- 
jected to external constraints (holonomic and nonholonomic) [3] . These systems 
have a wide application in many different areas: engineering, optimal control the- 
ory, mathematical economics (growth economic theory), subriemannian geometry, 
motion of microorganisms, etc. 

Constrained variational calculus have a rich geometric structure. Many of these 
systems usually exhibit invariance under the action of a Lie group of symmetries 
and they can be notably simplified using their symmetric properties by reducing 
the degrees of freedom of the original system. In previous studies it is imposed a 
separate study for each class of systems since the lack of a unified framework for 
dealing simultaneously with all the systems. 

Recent investigations have lead to a unifying geometric framework to covering 
these plethora of particular situations. It is precisely the underlying structure of a 
Lie algebroid on the phase space which allows a unified treatment. This idea was 
first introduced by Weinstein [34j in order to define a Lagrangian formalism which 
is general enough to account for the various types of systems. The geometry and 
dynamics on Lie algebroids have been extensively studied during the past years. In 
particular, in [27], E. Martinez developed a geometric formalism of mechanics on Lie 
algebroids similar to Klein's formalism of the ordinary Lagrangian mechanics and, 
more recently, a description of the Hamiltonian dynamics on a Lie algebroid was 
given in [531I2H]. The key concept in this theory is the prolongation, T^E, of the Lie 
algebroid over the fibred projection r (for the Lagrangian formalism) and the pro- 
longation, 7^E*, over the dual fibred projection t* : E* Q (for the Hamiltonian 
formalism). See [53] for more details. Of course, when the Lie algebroid is E = TQ 
we obtain that 7^ E = T{TQ) and J'^E* = T{T*Q), recovering the classical case. 
Another approach to the theory was discussed in 1?!] . The existence of symmetries 
in these systems makes interesting to generalize the Gotay-Nester-Hinds algorithm 
[13] to the case of Lie algebroids with a presymplectic section. These results are 
easily extended to the case of implicit differential equations on Lie algebroids. 

The second author and collaborators analyzed the case of nonholonomic me- 
chanics on Lie algebroids 8]. Now, we also pretend to study singular Lagrangian 
systems and vakonomic mechanics on Lie algebroids (obtained through the appli- 
cation of a constrained variational principle). 

The paper is organized as follows. In Section j^] we recall the notion of a Lie 
algebroid and several aspects related with it. In particular, we describe the pro- 
longation 7^E of a Lie algebroid E over the projection t : E ^ Q and how to 
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use this construction to develop Lagrangian mechanics on a Lie algebroid E with 
a Lagrangian function L : E —^ M., introducing several important objects, such as 
the Lagrangian energy El and the Cartan 2-section uj^. When the Lagrangian 
function is regular (that is, lul is nondegenerate) , we have existence and unique- 
ness of solutions for the Euler-Lagrange equations. However, if the Lagrangian is 
singular (or degenerate) we cannot guarantee these results. Motivated by this fact, 
in Section [3] we introduce a constraint algorithm for presymplectic Lie algebroids 
which generalizes the well-known Gotay-Nester-Hinds algorithm. In addition, we 
show that a Lie algebroid morphism which relates two presymplectic Lie algebroids 
induces a relation between the two associated constraint algorithms. 

In Section [H we apply the results of Section [3] to singular Lagrangian systems 
on Lie algebroids. More precisely, given a Lie algebroid t : E —* Q and a singular 
Lagrangian function L : i? — > K, we look for a solution X of the presymplectic 
system {7^E, ujl, <f ^El) which is a SODE along the final constraint submanifold. 
An example for an Atiyah algebroid illustrates our theory. 

In Section O we develop a geometric description of vakonomic mechanics on Lie 
algebroids. In this setting, given a Lie algebroid t : E ^ Q, we have a pair {L, M) 
where L is a Lagrangian function on E and Af C iJ is a constraint submanifold. 
In Section [5.11 we deduce the vakonomic equations using our constraint algorithm 
and study the particular case when it stops in the first step. In this situation, 
if the restriction of the presymplectic 2-section to the final constraint algebroid 
is symplectic, one can introduce the vakonomic bracket, which allows us to give 
the evolution of the observables. On the other hand, it is well know that classical 
vakonomic systems can be obtained from a constrained variational principle. This 
can also be done for vakonomic systems on Lie algebroids, as it is shown in Section 
15.21 In the particular case when we do not have constraints, our approach can be 
seen as the Skinner- Rusk formulation of Lagrangian Mechanics on Lie algebroids. 
This is explained in Section 15.31 where we also illustrate our results with several 
interesting examples: If E is the standard Lie algebroid ttq : TQ — > Q, then we 
recover some well-known results (see [H]) for vakonomic systems; in the case when 
E = g, a. real Lie algebra of finite dimension, we are able to model a certain class 
of Optimal Control problems on Lie groups (see [lHIH]); for the Atiyah algebroid, 
we analyze some problems related with reduction in subriemannian geometry by 
means of the non-holonomic connection 4 . Finally, we study Optimal Control on 
Lie algebroids as vakonomic systems and, as an illustration of our techniques, we 
find the equations of motion of the Plate-Ball system. 



2. Lie algebroids 

Let E' be a vector bundle of rank n over a manifold Q of dimension m and 
T : E ^ Q he the vector bundle projection. Denote by r{E) the C°°(Q)-module 
of sections oi t : E ^ Q. A Lie algebroid structure (|-,-],p) on £' is a Lie 
bracket |-, •] on the space T{E) and a bundle map p : E ^ TQ, called the anchor 
map, such that if we also denote by p : r{E) — > X{Q) the homomorphism of 
C°°((3)-modules induced by the anchor map, then 

lX,fY]^flX,Yj+piX)if)Y, 

for X,Y € r{E) and / e C°°{Q). The triple [E, {■, is called a Lie algebroid 
over Q (see |25j). 

If {E, |-, -l, p) is a Lie algebroid over Q, then the anchor map p : T{E) X(Q) 
is a homomorphism between the Lie algebras {r{E), |-, •]) and (X(Q), [•, •]). 
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Standard examples of Lie algebroids are real Lie algebras of finite dimension and 
the tangent bundle TQ of an arbitrary manifold Q. 

Another example of a Lie algebroid may be constructed as follows. Let tt : P ^ Q 
be a principal bundle with structural group G. Denote by $ : G x P ^ P the free 
action of G on P and by r$ : G x TP ^ TP the tangent action of G on TP. Then, 
one may consider the quotient vector bundle Tp\G : TP/G Q = P/G and the 
sections of this vector bundle may be identified with the vector fields on P which 
are invariant under the action Using that every G-invariant vector field on P is 
TT-projectable and the fact that the standard Lie bracket on vector fields is closed 
with respect to G-invariant vector fields, we can induce a Lie algebroid structure 
on TP/G. The resultant Lie algebroid is called the Atiyah (gauge) algebroid 
associated with the principal bundle tt : P — > Q (see [23l [25] ). 

If {E, |-, is a Lie algebroid, one may define the differential of E, : 

r{A''E*) T{a''+^E*), as follows 

k 

4=0 

+ Y.{-ly+'^l{lx,,x,lXo,...,x^,■■■XJ,■■■,Xk), 

i<j 

for fj. e r(A'''£'*) and X^, ...,Xk G T{E). It follows that (d^)^ = 0. Moreover, if 
X e r(£'), one may introduce, in a natural way, the Lie derivative with respect 
to X, as the operator ilf : r(A'=£'*) ^ T{a''E*) given by = ix od^ + d^ oix. 

Note that ii E ^ TQ and X e r{E) = X{Q) then d^'^ and Z^'^ are the usual 
differential and the usual Lie derivative with respect to X, respectively. 

If we take local coordinates (x*) on Q and a local basis {ba} of sections of E, 
then we have the corresponding local coordinates (a;*,?/"^) on E, where y'^(e) is 
the A-th coordinate of e £ £^ in the given basis. Such coordinates determine local 
functions p\, C^g on Q which contain the local information of the Lie algebroid 
structure and, accordingly, they are called the structure functions of the Lie 
algebroid. They are given by 

d 

P{e-A) = Pa-q^ and |eA,eB] = GabGc- 
These functions should satisfy the relations 



dPB . dp\ _ c 



D 



cyclic{A,B ,C) 

which are usually called the structure equations. 

If / G G°°(Q), we have that 

where {e^} is the dual basis of {e^}. On the other hand, if 6* G T{E*) and 9 — 6ce^ 
it follows that 



In particular. 
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On the other hand, if {E, |-, p) and {E' , |-, -Y^p') are Lie algebroids over Q 
and Q', respectively, then a morphism of vector bundles {F, /) from E to E' 



E 



Q 



F 



f 



E' 



Q' 



is a Lie algebroid morphism if 

for e riA^i?')*). (2.1) 

Note that (F, f)*(j)' is the section of the vector bundle t\^E* — > Q defined by 

((F, /)*(/.').(ai, . . . , afc) = 0'^(,)(F(ai), . . . , F(a,)), 
for X € Q and ai, . . . , g F^;. We remark that (|2.1|1 holds if and only if 

° /) - (F, fY{d^'g'). for g' G C-(g'), 

= (F,/)*(d^'a'), for a' e r((F')*). 

If {F, f) is a Lie algebroid morphism, / is an injective immersion and F\e^ '■ 
Ex — > E'jf^^^ is injective, for all x G Q, then {E,\-,-\et Re) is said to be a Lie 
subalgebroid of (£", |-, -Je', Pe')- 

If Q = Q' and f = id : Q ^ Q then, it is easy prove that the pair {F, id) is a Lie 
algebroid morphism if and only if 

FIX, Yj - IFX, FYf, p'{FX) - 

for x,r e r(£;). 

Let {E, |-, •], p) be a Lie algebroid over a manifold Q and F* be the dual bundle 
to E. Then, E* admits a linear Poisson structure {■,•}£•, that is, 

{•, ■}e' ■■ C°°{E*) X ^^{E*) C°°{E*) 

is a E-bilinear map, 

{F, G}e' = -{G, F}e' , (skew-symmetry), 

{FF',G}e' F{F',G}£;. +F'{F,G}£;., (the Leibniz rule), 
{F, {G, }£. + {G, {il, F}b. lij. 

{F, GIb-Ib- 0: (the Jacobi identity), 

for F,F',G,H e C°°{E*) and, in addition, 

P, P' linear functions on E* ^ {P, F'}^;. is a linear function on E*. 

If (x*) are local coordinates on an open subset [/ of Q, {ca} is a local basis of r(F) 
on U and F,G G C°°{E*) then the local expression of the Poisson bracket of F and 
G is 



{F,G}e^=p\ 



fdF dG dF dG\ dF dG 

^abPc- 



(2.2) 



\ dx'- dp A dp A dx^ J dp A dps ' 

where {x'^^pa) are the corresponding coordinates on E* (for more details, see [25]). 
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2.1. The prolongation of a Lie algebroid over a fibration. Let (i?, 

be a Lie algebroid of rank n over a manifold Q of dimension m and tt : P — > Q be 

a fibration, that is, a surjective submersion. 

We consider the subset T^P oi ExTP defined by T^P = |J T^P, where 

peP 

7^P - {{b, v) G X TpP I pib) = {Tp7r)iv)} 

and Ttt : TP TQ is the tangent map to tt. 
Denote by t'^ : I^P — > P the map given by 

T^ib,v)^Tp{v), 

for (&, v) G T^P, Tp : TP P being the canonical projection. Then, if m' is the 
dimension of P, one may prove that 

dim 7p P = n + ni — m. 

Thus, we conclude that T^P is a vector bundle over P of rank n + m' — m with 
vector bundle projection : 7^ P P. 

A section X oi t'^ : T^P — * P is said to be projectable if there exists a section 
X oi T : E ^ Q and a vector field U on P which is 7r-projectable to the vector field 
p{X) and such that X{p) = {X{tt{p)), U{p)), for all p G P. For such a projectable 
section X, we will use the following notation X = {X, U). It is easy to prove that 
one may choose a local basis of projectable sections of the space r(T^P). 

The vector bundle : 7^ P P admits a Lie algebroid structure (|-, •J'^, p'^). 
In fact, 

[(Xi,C/i),(X2,C/2)r = (IXi,X2l,[C/l,C/2]), p'^(Xi,C/i) =C/i. 

The Lie algebroid (T-^P, |-, -J" , p") is called the prolongation of E over tt or the 
E-tangent bundle to P. Note that if pri : J^P E is the canonical projection 
on the first factor, then the pair {pri,7r) is a morphism between the Lie algebroids 
(T^P, I-, -r, p'') and {E, |-, ■], p) (for more details, see [13]). 

Example 2.1. Let (P, be a Lie algebroid of rank n over a manifold Q 

of dimension m and t : P — > Q be the vector bundle projection. Consider the 
prolongation 7^ E oi E over r, 

T^P = {(e, v)eExTE\ p{e) = {Tt){v)}. 

7^E is a Lie algebroid over E of rank 2n with Lie algebroid structure (|-, ■Y,p'^). 

If (a;*) are local coordinates on an open subset U oi Q and {e^} is a basis of 
sections of the vector bundle t^^{U) U, then {Xa, "^a} is a basis of sections of 
the vector bundle iT'^)~^iT~^{U)) — » t^^(C/), where t"^ : T-^P ^ P is the vector 
bundle projection and 

d 

lAie) = {eA{T{e)),p\ — ^), 

for e G T^^(P). Here, p\ are the components of the anchor map with respect to 
the basis {ca} and {x^,y"^) are the local coordinates on E induced by the local 
coordinates (x') and the basis {e^}- Using the local basis {Xa,^a}, one may 
introduce, in a natural way, local coordinates {x'^,y^; z"^,v^) on T^P. On the 
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other hand, we have that 

for all A and B, C'^g being the structure functions of the Lie bracket |-, •] with 
respect to the basis {e^i}. 

The vector subbundle {7^E)^ of 7^E whose fibre at the point e e E is 

{7fEf - {(0,1;) eEx T,E/ {T,t){v) - 0} 

is called the vertical subbundle. Note that {7^E)^ is locally generated by the 
sections {V^}. 

Two canonical objects on T^E are the Euler section A and the vertical 
endomorphism S. A is the section of 7^E E locally defined by 

A^V^Va, (2.3) 

and S is the section of the vector bundle {7^ E) ® (7^ E)* E locally characterized 
by the following conditions 

SXa = Va, SVa = 0, for all A. (2.4) 

Finally, a section ^ of T^S E is said to be a second order differential 
equation (SODE) on E if S{S,) = A or, alternatively, pri{£,{e)) = e, for all e G E' 
(for more details, see [23j). 

Example 2.2. Let {E, |-, •], p) be a Lie algebroid of rank n over a manifold Q of 
dimension m and r* : E* — > Q be the vector bundle projection of the dual bundle 
E* to E. 

We consider the prolongation 7^E* of E over t*, 

7^E* = {(e', v)eExTE*\ p(e') = (Tt*)(i;)}. 

7^E* is a Lie algebroid over E* of rank 2n with Lie algebroid structure (|-, - J^ , ). 

If (a;*) are local coordinates on an open subset U of Q, {ba} is a basis of sections 
of the vector bundle t^^{U) — > U and {e'*} is the dual basis of {ba}, then {^a, '^^} 
is a basis of sections of the vector bundle (t^ )^^{{'''*)^^{U)) (t*)^^(J7), where 
t'^ : 7^E* E* is the vector bundle projection and 

^^(e*) = (eA(r*(e*)),p^^A^^J, 

y^{e*) = (0,^ ), 

for e* e (T*)^^(t/). Here, {x'^,pa) are the local coordinates on E* induced by the 
local coordinates (x*) and the basis {e^} of T{E*). Using the local basis {y^i, J*"^}, 
one may introduce, in a natural way, local coordinates {x^ ,pA] m^, qa) on 7^ E* . If 
u)* is a point of (r'^ )^^{{'''*)^^{U)), then (x^pa) are the coordinates of the point 



T 



(lu*) e (r*)-i(C/) and 



tu* = «^yA(T"* (^*)) + <7a5'-^(t^* (c^*)). 
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On the other hand, we have that 

for all A and B. Thus, if {^^ , Ta} is the dual basis of {^a, J''^}, then 
/ - Paq^,"^ + dpA 

for / G C°°(i;*). 

We may introduce a canonical section As of the vector bundle {7^E*)* E* 
as follows. If e* G i?* and (e, w) is a point of the fiber of 7^E* over e*, then 

XE{e*){e,v) = e*{e). 

Xe is called the Liouville section of (T^ £'*)*. 

Now, the canonical symplectic section 51^; is the nondegenerate closed 2- 
section defined by 

Qe ^ -tf"^' Xe- (2.5) 

In local coordinates, 

XEix\pA) = PA^^, 

nE = y^AVA + le%pc^^ A y^, 

(for more details, see {53] ) . 

2.2. Lagrangian and Hamiltonian mechanics on Lie algebroids. Given a 
Lagrangian function L G C°°{E) we define the Cartan 1-section 9l G T{{7^ E)*), 
the Cartan 2-section uol G T{A^{7^ E)*) and the Lagrangian energy El G 

C^{E) as 

9l = S*{cf''^L), ujL = -(f''^OL and El = L1''^L-L. 

If (x*, y"^) are local fibred coordinates on E, {p\, C^^) are the corresponding local 
structure functions on E and {Xa, "^a} is the corresponding local basis of sections 
of 7^E then 

From (1131), (EH), (ESI) and jUIl), it follows that 

isxt^L ^ -S*{ixuJL), lAUJL = -S*{(f ^El), (2.8) 

for X G r(T^i;). 

Now, a curve t — > c(t) on _E is a solution of the Euler- Lagrange equations for 



El = - L. (2.7) 



L if 



c is admissible (that is, p{c{t)) — m{t), where m = t o c) and 
*(c(t),c(t))^L(c(0) - d^^^-Bi(c(t)) = 0, for all t. 
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If c{t) = (a;*(t), y'^it j) then c is a solution of the Euler-Lagrange equations for L if 
and only if 

■ i , A d / dL \ dL : dL ^ . ^, 

^ =pV, Tjrh^ +7O-eSsy''-PA77-=0. (2.9) 



dtKdy^J dyC '^''dx^ 

Note that if E is the standard Lie algebroid TQ then the above equations are 
the classical Euler-Lagrange equations for L : TQ R. 

On the other hand, the Lagrangian function L is said to be regular if cj^ is a 
symplectic section, that is, if is regular at every point as a bilinear form. In 
such a case, there exists a unique solution verifying 

i^^LOL - d^^'^EL = . (2.10) 

In addition, using ()2.8p . it follows that is(_,^LOL = i^^L which implies that 
is a SODE section. Thus, the integral curves of (that is, the integral curves of 
the vector field P^{^l)) are solutions of the Euler-Lagrange equations for L. is 
called the Euler-Lagrange section associated with L. 

^From (12. 6p , we deduce that the Lagrangian L is regular if and only if the matrix 
(Wab) = i^ Q g j is regular. Moreover, the local expression of is 

where the functions satisfy the linear equations 

^'"^ r/^ + l^P^y^ + l^eS^y^ - = 0, for all A. 



dyBQyA-' dx'dy^'^''" dy(^ """^ 

Another possibility is when the matrix (Wab) = (^ g aq b ) '^^'^ regular. This 

type of lagrangians are called singular or degenerate lagrangians. In such a 
case, is not symplectic and Equation (|2.10p has no solution, in general, and even 
if it exists it will not be unique. In the next section, we will extend the classical 
Gotay-Nester-Hinds algorithm [13 for presymplectic systems on Lie algebroids, 
which in particular will be applied to the case of singular lagrangians in Section |4l 

For an arbitrary Lagrangian function L : E ^ M., we introduce the Legendre 
transformation associated with L as the smooth map legL : E ^ E* defined by 

for e, e' £ Ex, where z' <E TfE E^ x T^E satisfies 

pri(z') = e', 

pri : T^E — > E being the restriction to J^E of the canonical projection pri : 
E X TE E. The map legL : E ^ E* is well defined and its local expression is 

legL{x\y^) = {x\g^). (2.11) 

The Legendre transformation induces a Lie algebroid morphism 

7legL:7^E ^7^E* 

over legL '■ E —>■ E* given by 

{7legL){e,v) = {e,{TlegL){v)), 
where TlegL ■ TE TE* is the tangent map to legL : E ^ E*. 
We have that (see |23j ) 

{7legL,legLnXE) = eL, [7 legLMgL)* {^e) = o^l- (2.12) 
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On the other hand, from (|2.1ip . it follows that the Lagrangian function L is 
regular if and only if legL '■ E E* is a local diffeomorphism. 

Next, we will assume that L is hyperregular, that is, leg^ : E E* is a 
global diffeomorphism. Then, the pair (TZeg^, legL) is a Lie algebroid isomorphism. 
Moreover, we may consider the Hamiltonian function H : E* defined by 

H = Elo legl^ 

and the Hamiltonian section G T{T^E*) which is characterized by the con- 
dition 

The integral curves of the vector field {^h) on E* satisfy the Hamilton 
equations for H 

dx"- j dH dpA ( I dH dH \ 

for i e {1, . . . , m} and A e {1, . . . , n} (see [53]). 

In addition, the Euler-Lagrange section associated with L and the Hamilto- 
nian section S^h are (T/egi, ^e^^,) -related, that is. 

Thus, if 7 : / ^ is a solution of the Euler-Lagrange equations associated 
with L, then /i — legL o 7 : / — > E* is a solution of the Hamilton equations for H 
and, conversely, if /i :/—>_£* is a solution of the Hamilton equations for H then 
7 = leg~j^^ o /X is a solution of the Euler-Lagrange equations for L (for more details, 
see 1231). 



3. Constraint algorithm and reduction for presymplectic Lie 

algebroids 

3.1. Constraint algorithm for presymplectic Lie algebroids. Let r : — > Q 

be a Lie algebroid and suppose that Vt S T{/\^E*). Then, wc can define the vector 
bundle morphism \>q : E ^ E* (over the identity of Q) as follows 

t'o(e) — i{e)n{x), for e £ E^- 

Now, if a; G Q and is a subspace of E^ , we may introduce the vector subspace 
of Ex given by 

F^ ^{eeE^l n{x){e, /) = 0, V/ G FJ. 

Then, using a well-known result (see, for instance, |24j ) . we have that 

diuiF^ = dimEx - dimFj, + dim{E^ F^). (3.1) 

On the other hand, if bjij. = ^ii\E^ it is easy to prove that 

bo.(F.)C(F,^)0, (3.2) 

where (F^)° is the annihilator of the subspace F^ . Moreover, using (|3.ip . we 
obtain that 

dim{F^f = diniF^ - dim{E^ D F^) = dim{\>nSFx))- 
Thus, from p.2p . we deduce that 

\>nAF.)^{F,^r- (3.3) 
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Next, we will assume that £7 is a presymplectic 2-section (d^fi = 0) and that 
a e r{E*) is a closed 1-section {d^a — 0). Furthermore, we will assume that the 
kernel of SI is a vector subbundlc of E. 

The dynamics of the presymplectic system defined by (fi, a) is given by a section 
X G r(i5) satisfying the dynamical equation 

ix^^a . (3.4) 

In general, a section X satisfying p.4p cannot be found in all points of E. First, 
we look for the points where (|3.4p has sense. We define 

Qi ^ {x e Q\3e e E^ : i{e)n{x) = a{x)} = {x G Q | a{x) e {Ex)}- 

^From p.Sp . it follows that 

Qi^{x£Q\a{x){e) = 0, for all e £ Kern{x) = E^}. (3.5) 

If Qi is an embedded submanifold of Q, then we deduce that there exists X : 
Qi ^ E a section oi t : E Q along Qi such that (|3.4p holds. But p{X) is not, 
in general, tangent to Qi. Thus, we have that to restrict to Ei = p^^{TQi). We 
remark that, provided that Ei is a manifold and ti = t\Ei : -Ei — > Qi is a vector 
bundle, ti : Ei ^ Qi is a Lie subalgebroid oi E ^ Q. 

Now, we must consider the subset Q2 of Qi defined by 

= {xeQi\ a{x){e) = 0, for all e e (i;i);^ = (p-i(r,Qi))^}. 

If (52 is an embedded submanifold of Qi, then we deduce that there exists X : 
Q2 ^ El a. section of ti : £'1 ^ Qi along Q2 such that (|3.4p holds. However, 
/9(X) is not, in general, tangent to Q2- Therefore, we have that to restrict to 
E2 — p^^{TQ2)- As above, if T2 = t\e2 : £'2 — > Q2 is a vector bundle, it follows 
that T2 : £2 — * Q2 is a Lie subalgebroid of ti : Ei —* Qi. 

Consequently, if we repeat the process, we obtain a sequence of Lie subalgebroids 
(by assumption): 



fe+i ^ 
Ek+1 ^ 



Qk - 

Ek 



Q2 ' 

T2 

E2 ' 



Tl 

El ' 



To = T 
i?0 ~ -E 



where 



and 



Qk+i = [x £Qk\ a{x){e) 



Ek+i=p-'{TQ 



0, forallee(p-i(T,Qfe))^} 

-1/ 



fe+ij 



(3.6) 



If there exists k € N such that Qk = Qfc+i, then we say that the sequence 
stabilizes. In such a case, there exists a well-defined (but non necessarily unique) 
dynamics on the final constraint submanifold Qf = Qk- We write 



Qf — Qk+i = Qk, Ef = Ek+i — Ek — p ^(TQk). 

Then, Tf = Tk : Ef = Ek ^ Q f = Qk & Lie subalgebroid of r : — > Q (the Lie 
algebroid restriction oi E to Ef). From the construction of the constraint algorithm, 
we deduce that there exists a section X G r{Ef), verifying p.4p . Moreover, if 
X G ^{Ef) is a solution of the equation p.4p . then every arbitrary solution is of 
the form X' = X + Y, where Y G r(£:/) and Y{x) G kerf7(a;), for all x GQf. In 
addition, if we denote hy Qf and a/ the restriction of SI and a, respectively, to the 
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Lie algebroid Ef — > Qf, we have that 0/ is a presymplectic 2-section and then 
any X G ^{Ef) verifying Equation p.4p also satisfies 

ix^f = aj (3.7) 

but, in principle, there are solutions of (|3.7p which are not solutions of (|3.4p since 
kerJl riEf C kerJl/. 

Remark 3.1. Note that one can generalize the previous procedure to the general 
setting of implicit differential equations on a Lie algebroid. More precisely, let 
T : i? — > Q be a Lie algebroid and 5* C be a submanifold of E (not necessarily a 
vector subbundlc). Then, the corresponding sequence of submanifolds of E is 

S'o ~ S 

Si =5onp-i(Tr(S'o)) 



'fc+i 



In our case, Sh 



'HTQk) (equivalently, Qk = T{Sk)). 



3.2. Reduction of presymplectic Lie algebroids. Let {E, |-, ■], p) and (£", |-, •]', 
p') be two Lie algebroids over Q and Q', respectively. Suppose that G T{a'^E*) 
(respectively, 51' £ T{/\'^{E')*)) is a presymplectic 2-section on t : E ^ Q (respec- 
tively, t' : E' Q') and that a e r{E*) (respectively, a' e r{{E')*}) is a closed 
1-section on t : E Q (respectively, t' : E' ^ Q'). Then, we may consider the 
corresponding dynamical equations 

ix^ = a, XeTiE), 

ix'^' = a', X' e T{E'). 

If we apply our constraint algorithm to the first problem, we will obtain a se- 
quence of Lie subalgebroids oi t : E ^ Q 



Qk+i ^ 

Tk+l 

Ek+1 ^ 



Qk - 

Tk 

Eh ^ 



Q2 ' 
Eo ' 



Qi ^ Qo — 



El ' 



To 



En = E 



In a similar way, if we apply our constraint algorithm to the second problem, we 
will obtain a sequence of Lie subalgebroids of t' : E' ^ Q' 



Q'k+i ^ Q'k 



'k+l 



K+1 



E'. 



E'2 ^ E[ 



^E'o = E' 



On the other hand, it is clear that the restriction il^ (respectively, 51'^,) of fl 
(respectively, fl') to the Lie subalgebroid : Ek ^ Qk (respectively, t'^ : E'^. ^ Q'^.) 
is a presymplectic section of : Ek Qk (respectively, t'/. : E'^, ^ Q'^,). Moreover, 
if ak (respectively, aJc) the restriction of a (respectively, a') to Tk ■ Ek ^ Qk 
(respectively, r(. : E[. — > Q'k), we may consider the dynamical problem 

iXk^k = ak, Xk e T{Ek), 

(respectively, ix' ^'k = a'fc, Xj^ e r(£'^)) on Tk : Ek Qk (respectively, t^. : E',^ ^ 

Q'k)- 



CONSTRAINED MECHANICS ON LIE ALGEBROIDS 



13 



Now, suppose that the pah (11, tt) 

n 

E E' 

T t' 

TT 

Q " Q' 

is a dynamical Lie algebroid epimorphism between E and E'. This means 
that: 

(i) The pair (11, tt) is a Lie algebroid morphism, 

(ii) TT : Q ^ Q' is a surjective submersion and Hi^^ : E^ E'^^^^^ is a Unear 
epimorphism, for all x € Q, and 

(iii) (n, Tr)*n' = n and (H, 7r)*a' = a. 

Then, we will see that the Lie subalgebroids in the two above sequences are 
related by dynamical Lie algebroid epimorphisms. First, we will prove the result 
for fc = 1. 

Lemma 3.2. // (H, tt) is a dynamical Lie algebroid epimorphism, we have that: 

a) 7r(Qi) = Q[ and n{Ei) = E[. 

b) If Xi e Qi, then 7r^^(7r(xi)) C Qi and Ker{Il\E^J ^ {Ei)^^. 

c) // TTi : Qi ^ Q'l and Hi : Ei E[ are the restrictions to Qi and Ei 
of TT : Q Q' and H : E E' , respectively, then the pair (HijTTi) is a 
dynamical Lie algebroid epimorphism. 

Proof. If a; e Q then, using that (H, 7r)*il' = fl and the fact that II\e^ ■ Ex E'^{x) 
is a linear epimorphism, we deduce that 

Tl{KerVt{x)) = KerVl' [nix)). (3.8) 
Thus, from ([33]) . ([XS]) and since (H, it)* a' = a, it follows that 

7r(Qi) c g;. 

Conversely, if x\ e Q\ and x £ 7r^^(x']^) then, using again (|3.5p . p.Sp and the 
fact that (n, 7r)*a' = a, we obtain that x e Qi. This proves that 

Q'l c ^(Qi) 

and the following result 

xi e Oi ^ 7r-i(7r(a;i)) C Qi. (3.9) 
Now, we will see that tti : Qi ^ Q\ is a submersion. 

In fact, if x\ £ Q\ and x'j^ = 7ri(a;i) then there exist an open subset V of Q' , 
x'l G U' , and an smooth local section s' : J7' — > Q of the submersion tt : Q — > Q' 
such that s'{x[) = xi. Note that, using p.9p . we conclude that the restriction 
s[ of s' to the open subset U[ — U' H Q[ of Q[ takes values in Qi. Therefore, 
s'l '. U{ C Q'^ — > Qi is a smooth map, s'i{x'i) — xi and s'j^ o tti — Id. Consequently, 
T^i '■ Qi Q'l is a submersion. 

Next, we will prove that 

mEi)x,) = n(p-i(r,,Qi)) = {p')-\T^ix,)Q[) - fo'^ e Qi. 
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Since 

(p'on)ii,^^ =(r7rop)|s^^, (3.10) 

it follows that 

n((i?i).J C {E[U,,y 

Conversely, suppose that e'l e {E[)^(^^^-^ = {p')~^{Tj^(^^^-^Q[). Then, we can 
choose e € E^^ such that n(e) — e\. Thus, from (|3.10p . we have that 

(T7r)(p(e))er,(,,)g;. 

Now, using that tti : Qi — > Q'x is a submersion, we deduce that there exists 
I'l G Txj^Qi such that 

{T7r){v,) ^ (T7r)(p(e)), 

that is, 

vi - p(e) e Tj;i(7r"^(7r(2;i))) C T^^Qi. 
Therefore, p{e) G T^^Qi and 

eep-i(r,,Qi)-(^i),,. 

On the other hand, if e G Ker{]l\E^-^): then 

o = n(e) e 

and, proceeding as above, we conclude that e g 

Finally, using that the pair (H, tt) is a dynamical Lie algebroid epimorphism, we 
obtain that the pair (IIi, tti) is also a dynamical Lie algebroid epimorphism. 

□ 

Next, we will prove the following theorem. 

Theorem 3.3. Let (n,7r) he a dynamical Lie algebroid epimorphism between E 
and E' . Then, we have that: 

(i) 7r(Qfe) = Q'f^ and n{Ek) = E'^, for all k. 

(ii) Ifxk G Qk, then 7r~^(7r(a;fc)) C Q^. and Ker{Il\E^J ^ {Ek)x^, for all k. 

(iii) // TTfe ; Qk Q'^ and Ilk '■ Ek E'f. are the restrictions to Qk and Ek 
of t: : Q ~> Q' and II : E E' , respectively, then the pair {Ilk,TTk) is a 
dynamical Lie algebroid epimorphism, for all k. 

Proof. The result holds for fc — 0, 1. Then, we will proceed by induction. 

Assume that the result holds for k G {0,1,..., N}. Then, we will prove it for 
k = N + 1. 

Note that if fc e {0,1,..., N} and Xk G Qk then, using the following facts 

{n,TTyn' - n, mEkU) - and n{E,j = 

we obtain that 

mEk)i) = iE'k)^i.^y (3.11) 
Thus, proceeding as in the proof of Lemma we deduce the result. 

□ 

We remark that the behavior of the two constraint algorithms is the same. In 
fact, if we obtain a final Lie subalgebroid Tf = Tk : Ef — Ek Q f = Qk for the first 
problem (that is, if Qk = Qk+i) then, from Theorem l3.31 it follows that QJc = Q'k+i 
and we have a final Lie subalgebroid Tj = rj^ : Ej — E'^ ^ Q'^ = Q'^ for the second 
problem. Conversely, if the second constraint algorithm stops at a certain k (that 
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is, Q'j, = Q'l^^i) then, using (|3.6p and p. lip , we deduce that Qk = Qk+i, i-e., the 
first constraint algorithm also stops at the level k. 

Now, suppose that X : Qk ^ Ek is a section of the Lie algebroid : Ek Qk 
such that 

and X is (11^, 7rfc)-projectable, i.e., there exists X' G r(i?^) satisfying 

X' O TTfe = life o X. 

Then, using that (11, 7r)*f7' = Vl and that (H, 7r)*a' = a, we obtain that 

In others words, X' is a solution (along Q'j,) of the second dynamical problem. 
Conversely, if X' E r(-B^) is a solution of the dynamical equation 

then we can choose X G r(i?fe) such that 

X'oTTk^IlkoX 

and, since (11, 7r)*ri' — H. and (n,7r)*Q;' = a, we conclude that 

4. Singular Lagrangian systems on Lie algebroids 

Let L : i? ^ M be a Lagrangian function on a Lie algebroid t : E Q. 

Denote hy lol and El the Cartan 2-section and the Lagrangian energy, respec- 
tively, associated with L. Then, ojl is not, in general, a symplectic section and, 
thus, the dynamical equation 

ix^^L = a El 

has not, in general, solution. Moreover, if there exists a solution of the above 
equation, it is not, in general, a second order differential equation and it is not, in 
general, unique. 

Note that the Legendre transformation legL : E ^ E* associated with L is not, 
in general, a local diffeomorphism. 

Definition 4.1. The Lagrangian function L is said to be almost regular if the 

following conditions hold: 

i) The subset Afi = legL{E) of E* is an embedded submanifold of E* . 

ii) The map legi : E — > Mi induced by the Legendre transformation is a 
submersion with connected fibres. 

In what follows, we will assume that L is an almost regular Lagrangian. 
Then, we may prove the following result. 

Proposition 4.2. The Lagrangian energy El is a basic function with respect to 
the submersion legi : E —> Mi, that is, there exists a Hamiltonian function H on 
Ml such that 

H o legi = El- 
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Proof. Suppose that e is a point of E and that (a;*, y"^) are fibred local coordinates 
in an open subset of E which contains to e. Then, using ()2.1ip . we deduce that 
^ = -^'^ofa'i 6 ^e-E is vertical with respect to the submersion legi : E Mi if 
and only if 

^^TT^Tb =0' for alls. 
Thus, \i X is vertical, from (|2.7p . it follows that 

dL , . R_ d'^L 



x{El) = y 



dy^ 



dy^dy^ 



dy^ 



= 0. 



This ends the proof of the result. 



□ 



Now, since r*^^^^ : Mi Q is a, fibration, one may consider the prolongation 
7^ Ml of the Lie algebroid t : E ^ Q over r*^^^ or, in other words, the £^-tangent 
bundle to Mi. 

Then, the canonical symplectic section He ox\ 7^E* E* induces a presym- 
plectic section Hi on the Lie algebroid 7^ Mi Mi. Moreover, the submersion 
legi : E Mi induces, in a natural way, a Lie algebroid epimorphism 

7legi : 7^ E 7^ Mi 



(over legi). In addition, using (|2.12p and Proposition 14.21 we obtain that 
{7legi,leginni) = ljl, {7 legi, legi)* {d''" El) = (f'^'^H. 
Therefore, we have proved the following result. 

Proposition 4.3. The pair (7 legi, legi) is a dynamical Lie algebroid epimor- 
phism between the presymplectic Lie algebroids [7^ E, uj^, d7 ^ El) and {7^ Mi, fii, 
(f^'^^H). 

The following diagram illustrates the above situation. 

7leqi , , 
— T^Afi 




Now, we consider the following dynamical equations 



ixojL = cF ^El, with X e r(T^£;) 



and 



iyHi = d'^'^^^^H, with Y e r(T^Mi). 



(4.1) 



(4.2) 



In general, a section X e T{7^E) (respectively, Y e T{7^Mi)) satisfying (|4.ip 
(respectively, (14. 2p ) can not be found in all the points of E (respectively. Mi). 
Thus, we must apply the general constraint algorithm developed in Section [3. II for 
an arbitrary presymplectic system. 
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Assume that this algorithm stops at the level k for the first dynamical equation, 
that is, there exists a Lie subalgebroid {'I^E)k of T^E over a submanifold Ek of 
E and a section Xk G T{{7^E)k) such that 

Note that {7^E)k = {p'')-^{TEk), where : T^E' TE is the anchor map 
of the Lie algebroid 7^E E. Moreover, using Proposition l4.3l and the results of 
Section [XH we deduce that the constraint algorithm also stops at the level k for 
the second equation. In fact, we have that: 

(i) legi{Ek) — Mk+i is a submanifold of Mi and 

(T^Mi)fe = {7legi){{7^E)k) = (p-^'r^TAh+i) 

is a Lie algebroid over Mk+i, : T^E* TE* being the anchor map of 
the Lie algebroid 7^E* E*. 

(ii) Uek G Ek then leg^\legi{ek)) C Ek and Ker{7^Jegi) C {7f^E)k. Note 
that, from (fT^ and (HHI]), it follows that 

KeriJeJegi) = KerLOLick) H (Tf^S)^. (4.3) 

(iii) If legk+i ■ Ek Af^+i and Legk+i ■ {7^E)k {7^Mi)k are the restric- 
tions to Ek and {7^E)k of legi : E = Eq ^ Mi and 7legi : 7'^E 
7^Mi, respectively, then the pair {Legk+i,legk+i) is a dynamical Lie al- 
gebroid epimorphism. 

(iv) If Xk e T{{7^E)k) is such that {tx,u;L)\E, = {d'^^^EL)\E, and Xk is 
(iegfe+i, Zeg/j+i)-projectable, i.e., there exists Yfe G r{{7-^ Mi)k) satisfying 

Yfe o legk+i = Legk+i o Xk, 

then 

(v) If Yk E T{{7^ Mi)k) is a solution of the dynamical equation 

then we can choose Xk E r{{7^ E)k) such that 

Yfe o ^e^fe+i = Legk+i o Xfe and {ixk^L)\Ek = [d^ ^EL)\Ek- 
Now, suppose that X E T{{7^ E)k) is a solution of the dynamical equation 

and that X is (Legfe+i, ^egfe+i)-projectable over Y E r((T^Mi)fe). 
Then, 

(*f1^i)|m,+, = (rf^"'^^H)|A/.+,. (4.4) 
The following diagram illustrates the above situation. 
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If pI : {'7^E)k -> TEk is the anchor map of the Lie algebroid {7^E)k Ek 
then the integral curves of the vector field pl.{X) don't satisfy, in general, the Euler- 
Lagrange equations for L. The reason is that the section X is not, in general, a 
SODE along the submanifold Ek of E. In other words, X doesn't satisfy, in general, 
the equation 

SX = /\\E,. 

A solution for the above problem is given in the following theorems. 
Theorem 4.4. (i) The subset of Ek defined by 

5^ = {e e Ek/{SX){e) = A(e)} 

is a submanifold of Ek ■ 

(a) There exists a Lie subalgebroid of (7^ E)k Ek (over S-^ ) such that if 
Leg^x : — > (T^Mi)^ and leggx : S'^ — > Mk+i are the restrictions of Legk+i ■ 
{7^E)k {7-^Mi)k and legk+i : Ek Mk+i to and S-^ , respectively, then 
the pair {Legj^x , leggx) is a Lie algebroid isomorphism. 

Theorem 4.5. There is a unique section G r(A^) satisfying the following 
conditions 

ii^xu;L)\s- = (cP"^El\s-, S{e) = A|sx. 

Theorem 4.6. If Pa^ ■ ^ TS*^ is the anchor map of the Lie algebroid 

S'^ then the integral curves of the vector field PA^i^'^) on S'^ are solutions of the 

Euler- Lagrange equations for L. 

Proof of Theorem \4-4\ We consider the smooth map Wx ■ Ek —>■ E defined by 

Wxie)^pri{X{e)), 

where pri : 7^E E is the restriction to T^E of the canonical projection pri : 
E xTE ^ E on the first factor. 

Now, we will proceed in several steps. 

First step: We will prove that 

n leg-l^{legk+iie)) - {Wx{e)}, for all e £ Ek- (4.5) 
E ee Ek then, using (g^, and the fact that 

we deduce that 

{SX - A)(e) e KerujL{e)- 
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On the other hand, it is clear that 

(5X-A)(e)e(Tfi?)^. 

Thus, from (|4.3p . we obtain that 

(5X-A)(e) e KeriT-Jegi). 

Therefore, if : 7^E TE is the anchor map of the Lie algebroid T^E, then 

pmSX - A)(e)) e Ker{T,legi) = Ker{TJegk+i). 

This imphes that X* = p^{SX — Ai^^.) is a vector field on Ek which is vertical 
with respect to the submersion legk+i : Ek —>■ Mk+i- 

Suppose that the local expression of the section X is 

X = X'^Xa + V^Va- (4.6) 

Then, 

SX-/\ = {X^-y^)VA and X* = {X^ - y^)-^. (4.7) 

Moreover, since X is (Leg/c+i, ^e(7fc_|_i)-projectable, the functions X^ are constant 
on the fibres of legk+i- Consequently, if e = [x^.y^) € E^, it follows that the 
integral curve oi X* over the point e is 

s^cT{s)^{x\,,X^ + e-\y^-X^)). 

In particular, 

<7[s) = (xlX^ + e-\y^ - X^)) € leg^l^{legk+i{e)), for aU s G R, 
which implies that 

lim a{s) EE {xl,X^) e leg^l^{legk+i{e)). 

S — *CX} 

Now, from ()4.6|) . we have that 

Wxie)^pn{X{e)) = {xlX^). (4.8) 
In addition, using (14. 7p and (14. 8p . one deduces that 

X*{Wx{e)) = 0. 

On the other hand, if e' G 5"^ n leg^^^{legk+i{e)) then, it is clear that 

r(e)=r(e') 

and, thus, e' = (xq, ?/"^(e')). Furthermore, since e' G 5*^, we obtain that 
= X*ie') = {X^-y^ie'))j^ 

which implies that 

e' ^{xl,X^)^Wx{e). 
Second step: We will prove that 

where Wx ■ Mk+i Ek is a. section of the submersion legu+i ■ Ek Mk+i- 
Therefore, 5'''^ is a submanifold of Ek and dim S-^ = dim Alk+i- 

In fact, suppose that e, e' G Ek and 

legk+i{e) = legk+i{e'). 

Then, 

Legk+i{X{e)) = Yilegk+i{e)) = Y{legk+iie')) = Legk+i{X{e')). 
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Consequently, 

Wx{e) = pri{X{e)) = pri{X{e')) = Wx{e'). 

So, we have that there exists a smooth map Wx ■ Mk+i — > Ek such that the 
foUowing diagram is commutative 

Wx 




Wx 
Mk+i 

Moreover, using ()4.5|) . we deduce that Wx ■ Mk+i — > -Efe is a section of the 
submersion legk+i ■ Ek — > Mk+i and 

S"" ^Wx{Mk+i). 

Third step: We wiU prove the second part of the theorem. 

The section Wx '■ Mk+i Ek induces a map 

7Wx ■■ (T^Mi)fc = {p-')-\TMk+i) ^ (T^-B)fe = {pn-\TEk) 

in such a way that the pair (TWx, Wx) is a Lie algebroid monomorphism. We wiU 
denote by the image of (T^A/i)^ by the map 7Wx- Then, it is clear that A-^ is 
a Lie subalgebroid (over S^) and the pair {7Wx, Wx) is an isomorphism between 
the Lie algebroids A-^ S-^ and {7^Mi)k Mk+i- In fact, the inverse morphism 
is the pair [LegAx ,leggx). 

□ 

Proof of Theorem \4.5\ We consider the section G r(y4'^) defined by 

^7Wx oY olegsx. 

Using (|4.4p and the fact that 

LegAX o^"" =Yolegsx, (4.9) 

it follows that 

{i^xUJL)\SX = {d'^''^EL)\sx. 
Now, from (|4.9p . we have that 

(e^ - X){e) e Ker{Legk+i)\i^7BE), = Ker{7Jegi), for all e G 5^, 

which implies that (see (|4.3p ) 

(S'^^)(e) = (5X)(e) = A(e), for aU e e 5"^. 

Next, suppose that ry is another section of the Lie algebroid A-"^ ~> S-^ such that 
(z,,cjl)|5x {<f''^EL)\sx, Sri = A|sx. 

Then, it is clear that 

iv - ^^)(e) e iferu;L(e) n {7f E)"^ , for all e G S"^, 
and, using ()4.3p . we deduce that 

- ?^)(e) G Ker{7Jegi) = ifer(Le5fc+i)|(T_E£;)^ , for all e 6 S"-^. 

Thus, 

LegA-{{r]~e){e)) = Q 
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and, since Leg^x : — > (T^Mi)^ is a vector bundle isomorphism, we conclude 
that 

□ 

Proof of Theorem The section is a SODE along the submanifold S-^ . 
Therefore, from (|2.3p and (j2.4p . we have that the local expression of is 

e = y''lB\S-+i'''^B\S-- (4.10) 

On the other hand, using (j2.6|l . (12. 7p . ()4.10p and the fact that {i(^xLOL)\sx = 
it follows that 

Now, the local expression of the vector field Pa^H^) on S-^ is 



Consequently, the integral curves of p^x (^^ ) satisfy the following equations 



which are the Euler-Lagrange equations for L. 

□ 

Remark 4.7. If we apply the results obtained in this Section for the particular 
case when the Lie algebroid E is TQ, we recover the results proved in [TTl [T^] for 
standard singular Lagrangian systems. o 

Example 4.8. To illustrate the theory we will consider a variation of an example 
of singular lagrangian with symmetry. This example corresponds to a mechanical 
model of field theories due to Capri and Kobayashi (see [3[7]). 
Consider the lagrangian function 

L = 1-1712 {±1 + yl) + y2X2 - ±22/2 ~x\~yl~xl~ y^. 



The configuration space is Q = M'* with local coordinates {xi,yi,X2,y2)- Clearly, 
the lagrangian is singular; in fact, since 

ujl = m2dx2 A dx2 + "^2^2/2 A dy2 + 2dy2 A dx2, 

then 

( d d d d 
Ker UJL = span <^ -— , — , -— , — 
[dxi dyi dxi dyi 

The system is invariant by the 5^ action 

S^xQ Q 

(a, (a;i,yi,X2,y2)) ' — > (xi, yi, 2:2 cosa - j/2 sina, a::2 sina + 7/2 cos a) 

Note that the action of 5^ on the open subset Q = x (M^ — {(0, 0)}) is free and 
then we may consider the reduced space Q/S^ and the Atiyah algebroid tq\S^ : 
TQ/S^ Q/S^. Taking polar coordinates X2 = pcos9 and y2 = psinO, we have 
that the canonical projection n : Q ^ M ^ Q / is given by 

'^{xi,yi,p,d) = {xi,yi,p). 
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It is clear that the Atiyah algebroid is isomorphic to the vector bundle TM x IR — » 
M. 

On the other hand, a local basis of S^-invariant vector fields on Q is { , , ^ , 
^}. These vector fields induce a local basis of sections {ei, 62, 63, eo} of the Atiyah 
algebroid tq\S^ : TQ/S^ ^ M = Q/S\ Moreover, if (|-, is the Lie algebroid 
structure on tq\S'^ : TQ/S^ M = Q/S^, it follows that |ei, ej] = 0, for all i and 
j, and 

Now, if (xi, 2/1, p, ii, yi, /9, r) are the local coordinates on TQ/S^ induced by the 
local basis {ei, 62, 63, eo}, then the reduced lagrangian is 

I = \mi {{? + {prf) + p'r -x\-y\- . 

Thus, the Euler-Lagrangc equations for / are: 

mip — {m-iT + 2)pr + 2p = 0, 

Tdirf? -\- f? = constant, 
x\ = 0, 
Vx = 0. 

The local basis {ei, 62, 63, eo} of r{TQ/S^) induces a local basis 
{Xi, X2, X3, Xo, Vi, ^^2, V3, Vo} 
of r(T^'3/'^'(TQ/S'^)). The presymplectic 2-section cji is written as 
ui = m2X^ AV^ + m2p^X° A V° - p(m2r + 2)X^ A X°. 
The energy function is 

El = \m2 {p^ + {prf) +xl+yl+p' 

and 

^t-Q/^Vq/s^^^ = m2pV3 + p2^2rV° + {m^pr'' + 2p)X3 + 2xtX^ + 22/1X2. 

Thus kerw/ = {Xi,X2,Vi,V2} and the primary constraint submanifold Ei d E = 
Eq is determined by the vanishing of the constraints functions: xi = 0, yi = 0. 
Now ((p'"ol^')-i(T£i))-L = kerw,, and therefore Ef = Ei. 

Any solution X S {7^E)i of the dynamical equation 

{ix^i\^, = {cP^''"^^'=^"'e,)\e, (4.11) 

is of the form: 

X = pXs + rXo + fVi + gV2 H V3 Vq, 

m2 m2p 

where / and g are arbitrary functions on Ei . 

The Legendre transformation legi is in this particular case: 

legi{xi,yi,p;xi,yi,p,r) = {xi,yi, p;0,0,m2p,m2p'^r + p"^) 

Therefore, the submanifold Mi = legi{TQ / S^) of T*Q/S^ is defined the constraints 
Pxi = and = where we choose coordinates {xi,yi, p;pxi,Pyi,Pp,Pr) on 
T*Q/S\ 
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In the induced coordinates {xi,yi, p;pp,Pr) on Mi, the haniihonian h : Mi M 

is: 

h{xi,yi,p;pp,pr) = ^^Pp + 2T7T.^p2 ~ ^^)^ +^i+yi + P^- 

Applying the constraint algorithm to the presymplectic Lie algebroid [7^ Mi, fli, 
^T^Mi^-j deduce that the final constraint submanifold M2 is determined by 

M2 = {{xi,yi,p;pp,pr) e Ml I xi = 0,2/1 = 0} 

where there exists a well defined solution of 

Now, we return to the lagrangian picture and wc study the SODE problem. 
Observe first that all the solutions X of Equation (|4.1ip are {Leg2, ?e(72)-projectable. 
Now, if we additionally impose the condition SX = Ai^;^^, that is, 

pVa + rVo = iiVi + yiVa + pVg + rVo 

along El, we obtain that the submanifold S-^ C Ei is uniquely defined as 

5^ = {(0,0, p; 0,0, p,r) eTg/51}, 

and, therefore, the section is the SODE defined by: 

f:^-nr ^rr , (^2r + 2) - 2p 2p(m2r + 1) 

5. VAKONOMIC MECHANICS ON LlE ALGEBROIDS 

5.1. Vakonomic equations and vakonomic bracket. Let r : i? ^ Q be a 
Lie algebroid of rank n over a manifold Q of dimension m and i : ^ M be a 
Lagrangian function on E. Moreover, let M d E he an embedded submanifold of 
dimension n + to — to such that tjv/ — t\m : M — > Q is a surjective submersion. 

Now, suppose that e is a point of M, TM(e) = x E Q, that (x*) are local 
coordinates on an open subset U of Q, x e [/, and that {ca} is a local basis of 
T{E) on U. Denote by (x',?/^) the corresponding local coordinates for E on the 
open subset t^^{U). Assume that 

M n T~\U) = {{x\y^) e r-i(C/) | ^"{x\y^) = 0, a = 1, . . . , m} 

where are the local independent constraint functions for the submanifold M. 
The rank of the (to x (n + TO))-matrix 

V dx^ ' dy"^ / 

is maximun, that is, to. On the other hand, since tm : Af — > Q is a submersion, 
we deduce that there exists Vi G T^M such that {TTM){vi) — gfrj^, for all i G 
{1,...,to}. Thus, 

d ^ d 



which implies that 

"a^ , for a e {1, ... ,to| and I e |1, ... ,m|. 

Therefore, the rank of the matrix 



V [g/ Q=l,...,m;A=l,...,n 
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is m. We will suppose, without the loss of generality, that the (to x TO)-matrix 

\dy^ |g/ Q = l,...,rn;S = l,...,rn 

is regular. Then, we will use the following notation 

for l<A<n, l<a<rfi and fh + 1 < a < n. 

Now, using the implicit function theorem, we obtain that there exist an open 
subset V of T~^{U), an open subset W C ^m+n-m ^^j^^ smooth real functions 

: ^ M, a = 1, . . . ,m, 

such that 

MnV = {{x\ y^)eV\y''^ y"^), a = 1, . . . , to}. 

Consequently, (x*, y") are local coordinates on M. We will denote by L the restric- 
tion of L to A/. 

Now, we will develop a geometric description of vakonomic mechanics on Lie 
algebroids, naturally generalizing the previous results of the third author and colla- 
borators P . Moreover, for the case M = E we also generalize the formulation given 
by Skinner and Rusk [32l[33] for singular lagrangians to general Lie algebroids. 

Consider the Whitney sum of E* and E, E* © E, and the canonical projections 
pri : E* ® E — > E* and pr2 : E* ® E — > E. Now, let Wq be the submanifold 
Wo = pr'^^{M) — E* y.Q M and the restrictions tti = pr i^t^^ and 1^2 = pr2\Wo- 
Also denote by v : Wq — > Q the canonical projection. The following diagrams 
illustrate the situation 




E* 



W2 



E 




Next, we consider the prolongation of the Lie algebroid E over t* : E* ^ Q (res- 
pectively, v : Wo Q). We will denote this Lie algebroid by 7^E* (respectively, 
7^Wq). Moreover, we can prolong tti : Wq ^ E* to a. morphism of Lie algebroids 
Ttti : T^Wo 7^E* defined by Ttti = {Id^Tni). 

If (x'^jPa) are the local coordinates on E* associated with the local basis {ca} of 
T{E), then {x\pA,y°') are local coordinates for Wq and we may consider the local 
basis {yA,?^,^,,} of r(T^Wo) defined by 

■ d 

yA(e*,e) = {eA{x),PA- 



d 



0), 



,0), 



op A |e. 

V,(e*,e) = (0,0,/^ ), 

where (e*, e) £ Wq and i^(e*, e) = x. If (|-, •]", p") is the Lie algebroid structure on 
T^Wo, we have that 

and the rest of the fundamental Lie brackets are zero. Moreover, 
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The Pontryagin Hamiltonian is a function in Wq = E* Xq M given by 

Hw„{e*,e) = (e*,e) - L(e), 

or, in local coordinates, 

Hwoi^^PA^y") =Pay'' +Pa^"{x\y'') - L{x\y''). (5.1) 

Moreover, one can consider the presymplectic 2-section VLq = {7iti^tti)*^e, 
where fi^ is the canonical symplectic section on 7^E* defined in Equation (|2.5p . 
In local coordinates, 

f^o ATA + ^eSsPcy^ Ay^. (5.2) 

Therefore, we have the triple {7^Wq, VLq, (P'^'^°Hwo) as a presymplectic hamil- 
tonian system. 

Definition 5.1. The vakonomic problem on Lie algehroids is find the solutions for 
the equation 

ixno = d''^^«Hwo, (5.3) 
that is, to solve the constraint algorithm for (T^T4^0i ^Oi ^"Hwo)- 

In local coordinates, we have that 

If we apply the constraint algorithm, 

Wi^ {w e E* Xq M\ d'^''^°Hw„HiY) = 0, Vr e Kerno{w)}. 
Since Kerflo — span{Va}, we get that Wi is locally characterized by the equations 

^a^d-'-^OH^M=Pa+pJ^~l^^-0, 



or 



dL 9*" 



Let us also look for the expression of X satisfying Eq. (|5.3p . A direct computation 
shows that 



X = J/'^ya + *"ya 

Therefore, the vakonomic equations are 



Pa 



Remark 5.2. We note that the vakonomic equations can be obtained following 
a constrained variational principle. Below, we will show this variational way of 
obtaining these equations (see Subsection l5.2p . o 
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Of course, we know that there exist sections X of 7^Wq along Wi satisfying 
(f531) . but they may not be sections of {p'')^'^{TWi) = T^VFi, in general. Then, 
following the procedure detailed in Section [3l we obtain a sequence of embedded 
submanifolds 

• • ■ ^ Wk+i --^ Wk ^ ■ . ■ ^ W2 ^ Wi Wa ^ E* xq M. 

If the algorithm stabilizes, then we find a final constraint submanifold Wf on which 
at least a section X G T{7^Wf) verifies 

As in [U], we analyze the case when Wf — W\. Consider the restriction ili of JIq 
toT^Wi. 

Proposition 5.3. fJi is a symplectic section of the Lie algehroid T^Wi if and only 
if for any system of coordinates (x"^ ,pA,y°') on Wq we have that 

drd^ ^^"drd^l ^ ^"'""^ 

Proof. It is clear that d^'^'^if^i = 0. 

On the other hand, if it; e Wi then, since the elements (c?^ ^"Lpa){'w) are inde- 
pendent in (T^Wo)*, we have that 

(note that dim (T^W"i)° = n — ffi). Moreover, using a well-known result (see, for 
instance, ^^), we deduce that 

dim (T^I^i)^'"" = dim {7^Wo)-dim (Tf I^i) -t-dim {KerQaHnl^i^Wi). (5.4) 

Now, suppose that fii is a nondegenerate 2-section on 7^Wi W\. Then, it is 
clear that 

Ker^^{w)r\7lW^ = {0}. 
Thus, the matrix (c?^ ^°i^a)(w)^h^w)) is regular, that is. 



\w 



dyaQyb " dy°-dy^ ■ 



Conversely, assume that the matrix (d^ Lpa){w){Vi,{'w)) is regular, then 

KernQ{w)r\7^,Wi^{Q}. 

This imphes that (see 

dim (Tf = dim {7^Wa) - dim (if I^i) = dim {Kerno{w)). 

Therefore, since iferf^o(w) C {7^Wi)-^'^^'> , it follows that 

Kernoiw) = (Tf 

and, consequently, 

(7fM^i)n(Tfl^i)^^"»={0}, 
that is, ili{w) is a nondegenerate 2-section. □ 

In what follows, we will use the following notation 

dL 92^° ^ „ 
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Remark 5.4. Suppose that the submanifold M is a real vector subbundle D (over 
Q) of E (note that t^j = r|£i : D — > Af is a surjective submersion). Then, we may 
consider local coordinates (x*) on an open subset U oi Q and a local basis {cq, Ca} 
of r{E) on U such that {ca} is a local basis of r(Z3). Thus, if (a;',y",i/°) are the 
corresponding local coordinates on t~^{U), we have that 

T^\U) = DnT-\U) EE {ix\y",y^) G T-\U)/y" = 0}. 

In other words, the local function is the zero function, for all a. Therefore, in 
this case, 

^ab ^ QyaQyb ' for aU fl and 6. 
We remark that the condition 

implies that the corresponding nonholonomic problem determined by the pair (i, D) 
has a unique solution (see [S]). o 

Remark 5.5. We remark that the condition det (3?ab) / implies that the matrix 
dipa \ 

— — ^ I is regular. Thus, using the implicit theorem function, we de- 

" / a.fc— m-f-l,...,n 

duce that {x^.pA^y"') are local coordinates for Wq on an open subset Aq C Wq in 
such a way that there exist an open subset W C R™+" and smooth real functions 



/i" : 14^ M, a = m + 1, . . . ,71, 

such that 

WinAo = {{x\pA,y'') G Ao/y" = /i°(x\pA), a = m + 1, . . . , n}. (5.5) 
Therefore, a local basis of r(T^W^i) is given by 

{y^l = (liA + PA%'^a)\W,.n = (J'-' + ^^^^-Mw.}- 

This implies that 

{yAl,a'^(Va)|H^J 

is a local basis of r(T^^Wo). Note that, from (|5.5p . it follows that (a;*,p^) are 
local coordinates on Wi. Moreover, if vi : Wi Q is the canonical projection and 
([•7 'FS P"^) is the Lie algebroid structure on 7^Wi Wi, we have that 

and the rest of the fundamental Lie brackets are zero. In addition, 

p-(y.i) = P^|„ p-(T-) = ^. (5.6) 

o 

Now, we will prove the following result. 

Theorem 5.6. If ili is a symplectic 2-section on the Lie algebroid 7^Wi — > Wi 
then there exists a unique section of T^Wi W\ whose integral curves are 
solutions of the vakonomic equations for the system {L,M). In fact, if Hy/i is the 
restriction to Wi of the Pontryagin Hamiltonian IIw„, then ^1 is the Hamiltonian 
section of Hwi with respect to the symplectic section ili , that is, 
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Proof. We have that (see (|5.ip ) 

Therefore, from (|5.2p . it follows that 
This implies that 

and, consequently, the integral curves of are solutions of the vakonomic equations 
for the system {L, M). 

Moreover, if e r{7^Wi) is another solution of the equation 
then it is also a solution of the equation 

which implies that Ci = ^ 



Suppose that {x'^,pa) are local coordinates on Wi as in Remark 15.51 and that 
{^^1,^5*} is the corresponding local basis of T(T''^Wi). Then, if {^f , Tai} is the 
dual basis of {^yii, '^t}: we have that (see 



f^i -yf ATAi + ie^cy^ Ayf. (5.7) 



On the other hand, from (15.61). it follows that 



Therefore, 



"Sai - {^abPc^ y PA—rr—)-'^! ■ 



dpA dpB dx 

Note that 

dx^ -\}dx^^ dx^dy'^^^'''^'^^)^^^' dpA -ydpA^ dpAdt^^^'''^^^,^^^ 
which implies that 

£,i[x^,Pb) - /i°(a:^^Ps)yai + *"(a;^Ai"(2;^PB))yai 

Now, we will introduce the following definition. 

Definition 5.7. The vakonomic system (L,M) on the Lie algebroid t : E ^ Q is 
said to be regular iffli is a symplectic 2-section of the Lie algebroid I^Wi — > Wi. 

Suppose that {L,M) is a regular vakonomic system and that Fi € C°°(W^i). 
Then, the Hamiltonian section of Fi with respect to fii is the section 'K^'^^ of 
the Lie algebroid 7^Wi W\ which is characterized by the equation 

i(5{^i)f^i=rf^"'^^^i- 

Note that 6 = "K^^ . 
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Using the Hamiltonian sections, one may introduce a bracket of functions {•, •}(l,m) • 
C°°(Wi) X C°°{Wi) ^ C°°{Wi) a.s follows 

for Fi,Gi e C°°{Wi). The bracket {■,-}{l,m) is called the vakonomic bracket 
associated with the system {L,M). 

Theorem 5.8. The vakonomic bracket {•, •}(L,Af) associated with a regular vako- 
nomic system is a Poisson bracket on Wi. Moreover, if Fi G C°°{Wi) then the 
temporal evolution of Fi, Fi, is given by 

Fi = {Fi, Hwt}(L,M}- 

Proof. Let b^i : T^Wi — > (T^Wi)* be the musical isomorphism induced by fii 
which is defined by 

boi(Xi) = forXiGTf^l^i andwieWi. 

Then, we may introduce the section Ai of the vector bundle A^T^Wi W\ as 
follows 

Ai(ai,/3i) = r!i(b^i oai,bs,i o/3i), for ai,/3i £ r(T^W^i)*. 

If [•, is the Schouten-Nijenhuis bracket associated with the Lie algebroid 7^W\ 
— > W\ then, using that d^^^^ili = 0, one may prove that 

IAi,Air^ =0. 

Thus, Ai is a triangular matrix, the pair (T^Wi, (T^PFi)*) is a triangular Lie 
bialgebroid in the sense of Mackenzie and Xu (see Section 4 in [H]) and {•, -{(l.m) 
is a Poisson bracket on W\ (see Proposition 3.6 in 

On the other hand, if F\ € C°°{Wi) and p'^^ is the anchor map of the Lie 
algebroid 7^Wi ^ Wi then 



Therefore 



Fi = {Fi, Hwi}(L,M)- 

□ 

Suppose that {x^,pa) are local coordinates on Wi as in Remark [5.51 and that 
{'^Ai,'?^} is the corresponding local basis of r(T^VFi). If i^i G C°^{Wi) then, from 
dSll), it follows that 

Consequently, using (|5.7p . we deduce that 

ru-ni dFi (J dF\_ , dFi . 

This implies that 



dFi dGi dFi dGi \ dFi dCi 



m,Gi}(i,M)-P^^^-^^j-e^5Pc^^. (5.8) 

Now, we consider the linear Poisson structure 

{•, ■}e' ■■ C°°{E*) X C°^{E*) C°°{E*) 

on E* induced by the Lie algebroid structure on E and the canonical symplectic 
structure fi^ on 7^ E* (see Section [2]). 
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Corollary 5.9. // (i, M) is a regular vakonomic system on a Lie algebroid E then 
the restriction {t:i)\Wi ■ ~^ ^* '"'i ■ ~* E* to Wi is a local Poisson 
isomorphism. Moreover, if 7{'Ki)\\y-^^ : 7^Wi I^E* is the corresponding prolon- 
gation then the pair ('J'(7ri)|n/j , (7ri)|vi/i ) local symplectomorphism between the 
symplectic Lie algebroids (T^VFijfii) and {7^ E* j^Ie). 

Proof. From Remark 15.51 we obtain that (tti)\Wi '■ ^ E* is a local difFeomor- 
phism. Furthermore, using (|2.2p and (|5.8p . we deduce that 

{L^ ° {t^i)\Wi,G ° {t^i)\Wi}{l,m) = {F,G}e- o {tti)\Wi, 

for F, G G C°° [E* ) . This proves the first part of the corollary. The second part 
follows from the first one and using the fact that (^'(Tri)^^^^ , {tii)\^^^)*Q,e = ili. □ 

Remark 5.10. If (7ri)ny^ : Wi £'* is a global diffeomorphism then we may 
consider the real function H on E* defined by 

In addition, if G T{7^E*) is the Hamiltonian section of H with respect to the 
symplectic section flE then and are {7{'Ki)\\y^, (7ri)|v^'J-related, that is, 

o 

5.2. The variational point of view. As it is well known, the dynamics in a 
vakonomic system is obtained through the application of a constrained variational 
principle [1]. More precisely, let Q be the configuration manifold, L : TQ K 
be the Lagrangian function and M C TQ be the constraint submanifold. Then, 
one can consider the set of twice differentiable curves which connect two arbitrary 
points a:, y G Q as 

G^{x, y) — {a : [to,ti] ^ Q | a is C^, a{to) = x and a{ti) — y}. 

This is a smooth manifold of infinite dimension and, given a G C^(x, y), the tangent 
space of C'^{x,y) is described as 

TaG^x, y)^{X: [to, ii] ^ TQ \ X is C\ X{t) G r,(t)Q, X{to) = and - 0}. 

The elements of TaG^{x,y) can be seen as the infinitesimal variations of the curve 
a. Next, one can introduce the action functional 6S : C^{x,y) R defined by 
a ^ SS{a) — J^^ L{a{t))dt. Thus, the associated vakonomic problem {Q,L,M) 

consists of extremizing the functional 5S restricted to the subset &{x, y), which is 
given by 

e{x,y) = {ae e\x,y)\a{t) G M„(t) = M D TQ\a{t)) , V t e [to,ti]}, 

that is, a G C^(a;, y) is a solution of the vakonomic system if a is a critical point of 
SS^Q2(^^ yY In the particular case when we do not have constraints, i.e. M = TQ, 
then we recover the variational way to obtain Euler-Lagrange equations. 

Now, consider a Lie algebroid t : E ^ Q and L : — > R a Lagrangian function 
on it. In [30j it is shown how to obtain the Euler-Lagrange equations (on the Lie 
algebroid) from a variational point of view. Let us recall some aspects related with 
this formulation. First, the set of i?-paths is defined by 

Adm{[to,ti],E) = {a : [<o,ii] E\poa= -^(r o a)}. 



CONSTRAINED MECHANICS ON LIE ALGEBROIDS 



31 



This set is a Banach submanifold of the set of paths in E whose base path is 
C^. Two iJ-paths ao and ai are said to be i?-homotopic if there exists a morphism 
of Lie algebroids 

^ :TI xTJ E, 

where I — [0, 1], J — [to,ti], a(s, t) = <&(<9t |(s.t)) and 6(s, t) ~ ^{ds\(s,t)), such that 

a(0,t) = ao(i), 5(s,io)-0, 
a(l,t) = ai(i), 5(s,ii) = 0. 

iJ-homotopy classes induce a second differentiable Banach manifold structure on 
Adm{[to,ti], E). The set of _E-admissible paths with this second manifold structure 
will be denoted by J'([io, ii], In addition, at each i?-admissible curve a, the 
tangent space is given in terms of the so-called complete lifts of sections. In fact, 

Tay{[toM.E)^{if eTaAdm{[toM].E)\T^{to)^Q and 77(^1) = 0}. 

(for more details, see [51 [30])- We recall that if {e^} is a local basis for E and ry is 
a time-dependent section locally given by 

then rf , the complete lift of 77, is the vector field on E given by 
V -V + (pb^t- - V ^AB)y 



With the second manifold structure that it is introduced on the space of i?-paths, 
it is possible to formulate the variational principle in a standard way. Let us fix 
two points x,y € M and consider the set J'([to, ti], of i?-paths with fixed base 
endpoints equal to x and y, that is, 

T{[to,h],E)y^{aen[to,ti],E)\T{a{tQ))=x and r(a(ti))=y}. 

In this setting, for the action functional SS : T'([io, ti], E) M. given by 

SS{a) = f ' L{a{t))dt, 



the critical points of SS on Tdtojti], E)^ are the curves a E J'([to, ^i], -B)^ which 
satisfy Euler-Lagrange equations (|2.9p (see [50]). 

Now, let (L, M) be a vakonomic system on the Lie algebroid t : E ^ Q. We will 
denote by 7{[to,ti], M)^ the set of -E-paths in M with fixed base endpoints equal 
to x and y 

n[to,ti],M)y = {aey{[to,h],E)y\a{t)eM,yte[to,h]}. 

We are going to consider infinitesimal variations (that is, complete lifts 77*^) tangent 
to the constraint submanifold M. We are going to assume that there exist enough 
infinitesimal variations (that is, we are studying the so-called normal solutions of 
the vakonomic problem). Since M is locally given by y" — '^°'{x'^,y"-) — 0, we 
deduce that the allowed infinitesimal variations must satisfy 

77'=(y"-^"(a;\y'')) = 0, 77(^0) = 0, 77(ti) = 0. 

Note that if a G ^([to, ti], M)^ then 

rj'^iy'^ -^"{x\y''))oa^O 

if and only if 

= p^77^_— + — ! — ^.Q'L yi^ri^ e^Ay^rj"^. 5.9) 

dt '^'^ ' dx' dt dy'' ""^^ ' 9y« ' ^ ' 
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Let US look for the critical points of the action functional SS : y{[to, ti], E)^ 
SS: Ti[to,hlE)y ^ M 



ait) 



If we consider our infinitesimal variations then, 



ds \s=o 



L{a,{t))dt 



to 



dL 

dL 

dx- 



dL 

dy 
dL 
dy 



L{a{t))dt. 



dL 9*" 



-Vn + t; TT—Vi 



dy"' dx' 



dl^d^ 



C>L ^ a. dL , a dL ^ 



Let Pa be the solution of the differential equations 



/ dL 

\dx^ '^^ dx^ 



where 



dL a*" 



dy"- ^" dy°- 



(5.10) 



Using Eq. ()5.9p . we get that 
d 



dt 



PaPaV 
dL 

dx 



dx^ ' ^" dt dy"^ ' ^"""^^ ' Qya 



Pc^'hAV^'v^ 



9*" 



If we use this equality, we deduce that 



ds |s=o 



L{a,{t))dt = 



Kdx'- 



9^ + d?(^"9^)'^ -^"^-^ 



9*" 



9y 



9y'' 



Finally, using Eq. (|5.10p and the fact that 
dr]'^ 



we obtain that 



ds |s=o 



to 



L{as{t))dt = 

/ di d^" 

\dx^ 9x* 



d / dL 



Pa 



dtKdy'^ dy" 
Since the variations if are free, we conclude that the equations are 



-'a-i 



Pa 



/ dl 9*'^ 
V9a;* ^'^ 9a::* 



)p\-V''^laPB 



9i 



9*" 



i 9?/° ^" 9y« 



V9x' dx"^ 



F^a-y'^abPE-^'^eLpB, 
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with Pa = — Pa^r, that is, we obtain the vakonomic equations for the vako- 
nomic system (L, M) on the Lie algebroid t : E ^ Q. 

5.3. Examples. 

5.3.1. Skinner- Rusk formalism on Lie algebroids. Suppose that t : E ^ Q is a, Lie 
algebroid, L : E ^M.\s a. Lagrangian function and M = E, that is, we do not have 
constraints. Then, 

Wii^E* ®E 

and the Pontryagin Haniiltonian Hwo '■ E* ® E ^M.is locally given by 
Hw„{x\PA,y^) ^ y'^PA ~L{x\y^). 

Let us apply the constraint algorithm. First, if fig is the presymplectic section on 
Wo given by — i!^WiTWi)*^E then, since Xerilo — span{VA\, the primary 
constraint submanifold Wi is locally characterized by 

VA = ^, (5.11) 

which are the definition of the momenta. W\ is isomorphic to E and the vakonomic 
equations reduce to 

- y^p\. 

Using (|5.1ip . it follows that Equations (|5.12p are just the Euler-Lagrange equations 
for L (see (|2.9p '). We remark that \i E — TQ, this procedure is the Skinner-Rusk 
formulation of Lagrangian Mechanics (see [321 [33]) ■ On the other hand, using 
Proposition 2131 we deduce that the vakonomic system is regular if and only if the 
Lagrangian function L is regular. 

5.3.2. The tangent bundle. Let E be the standard Lie algebroid ttq '■ TQ — > Q, 
M C TQ be the constraint submanifold such that ttq\m : M ^ Q is a surjective 
submersion and L : TQ ^ M be a standard Lagrangian function. Suppose that 
{q^, q^) = {q^, (t ^(f) local fibred coordinates on TQ and that the submanifold 
M is locally described by equations of the form 

= ^'"(g^,g''). 

As we know, the local structure functions of the Lie algebroid ttq '■ TQ Q with 
respect to the local coordinates {q^,q^) are 

pf = ,5f and Q% = Q. 

Thus, if we apply the results of Sections l5.1l and [5T^ to this particular case we recover 
the geometric formulation of vakonomic mechanics developed in [3]. In particular, 
the vakonomic equations reduce to 

f 9" = 

d / dL\ di _ . (9*° 
JtVd^) ' dq^ ^^"~d^ 
dL d-^^ 



r d /a^-^x 




Idt \ dq" ) 


dq" \ 
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5.3.3. Lie algebras. Let g be a real Lie algebra of dimension n. Then, g is a Lie 
algebroid over a single point. 

Now, suppose that £ is an afSne subspace of g modelled over the vector space 
C of dimension n — in and that eo € £, eo 7^ 0. We consider a basis {ba} = 
{ca, eo, eg} = {cq, Bq} of such that {cq} is a basis of C and 

[eA,eB] = G'^B^c- 

Denote by (i/°,y",y") = the linear coordinates on g induced by the basis 

{ca, Co, 65} = {cq, Cq}. Then, £ is given by the equations 

Next, assume that L : g — > R is a Lagrangian function and denote by L : £ — > R 
the restriction of L to £. Then, a curve 

a:t^ {y'^{t),y\t),y%t)) = (z/(t), 1, 0, ... ,0) 

in £ is a solution of the vakonomic equations for the constrained system (L, £) if 
and only if 



BT (5.13) 

Now, we consider the curve 7 in q* whose components with respect to the dual 
basis {e^} are 

oy |<j(t) 

that is, 7 is the sum of the curve ^ and the curve 

t ^ X{t) e C° 

whose components are X{t) = {0,pa{t)). Here, C° C g* is the annihilator of the 
subspace C. Then, a direct computation, using (j5.13p . proves that 7 satisfies the 
Euler-Poincare equations 

This result is just the "Optimization Theorem for Nonholonomic Systems 
on Lie groups^'' which was proved in 21j (see Theorem 5.1 in ^21j). 

5.3.4. Atiyah algebroids and reduction in subriemannian geometry. Let tt : P — > Q 
be a principal bundle with structural group G. The dimension of P (respectively, 
Q) is n (respectively, m). 

Suppose that 13 is a distribution on P such that 

TpP = Dp + VpH, for all p G P, 

where Vtt is the vertical bundle to the principal bundle projection tt. 

We also suppose that D is equipped with a bundle metric ( , )d. Assuming 
that both D and { , )d are G-invariant, then one may construct a nonholonomic 
connection as follows (see [1]). 

We will suppose that the space 

Sp = Dp n VpTT 

has constant dimension r, for all p E P. Under this condition, the horizontal space 
of the nonholonomic connection at the point p is 5"^^, where Sp is the orthogonal 
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complement of Sp on D with respect to the bundle metric { , )d- We will denote 
by w"'* : TP — *• g the corresponding Lie algebra-valued 1-form. 

Define a Lagrangian function on D hy L{vq) = ^{vq,Vq)D where Vq G Dg. 

Now, we consider the Atiyah algcbroid Tp\G : TP/G Q = P/G associated 
with the principal bundle n : P ^ Q = P/G. Note that the space of orbits D/G 
of the action of G on Z) is a vector subbundle (over Q) of the Atiyah algebroid. 

Next, we will obtain a local basis of T{TP/G) adapted to the vector subbundle 
D/G. 

For this purpose, we choose a local trivialization of the principal bundle tt : P ^ 
Q = P/G to he U X G, where U is an open subset of Q. Let e be the identity 
element of G and assume that there arc local coordinates (x*) on U. If i-^)^ is 
the horizontal lift of the vector field on U then (^r)'' is a G-invariant vector 
field and {{■£^)'^}i=i,...,m is a local basis of S-'-. 

On the other hand, we consider a family of smooth maps 

eA'-U^g, A = l,...,n-m, 

such that, for every q € U, {eA{q)} is a basis of g and {ea{q)}a=i,...,r is a basis of 
the vector space 

0('''^)={^G0|^p(g,e)GZ?(,,e)}. 
Here, is the infinitesimal generator of the free action of G on P associated with 

Now, we introduce the vertical vector fields on f/ x G given by 

Va : UxG T{U x G) TU x TG 

< 

{q,g) ^ eA{q){9) 

< 

eA{q) is the left-invariant vector field on G induced by the element eA{q) of g. 
Then, it is clear that 

is a local basis of G-invariant vector fields on P and 

is a local basis of the space of sections of the vector subbundle D ^ P. Thus, we 
have the corresponding local basis of sections 

of the Atiyah algebroid tp|G : TP/G Q = P/G which is adapted to the vector 
subbundle D/G^Q = P/G. 

Denote by (a;%i:%y",j/") the local coordinates on TP/G induced by the basis 
{ci, Ca, Ca}. Then D/G is locally characterized by the equations 

y" = o. 

On the other hand, if 

u^-'{^)=rt{q)eAiq), 

-Q§:^^ = Xwi<i)^c{q), [eA{q),eB{q)] = &ABec{q), 

for all g e J7, it follows that 

[ci, Bj] = "Bfjec, [e,, es] = /x-^ec, [ba, bb] = ^IbBc, 
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where 

arc dV^' 

■"ij ~ Q^j Qrj.i ^ » ^ J ^AB ^ -L t XjA -L J AiA: 

Moreover, if p : TP/G TQ is the anchor map of the Atiyah algebroid Tp\G 
TP/G -^Q = P/G, we have that 

d 



p{eA) = 0. 



Since the bundle metric ( , )d G-invariant, it induces a bundle metric ( , 

on D/G, with associated lagrangian function I : D/G R. Therefore, we deduce 

that a curve 

a:t^ {x\t),xHt), y\t),y^{t)) = (x^t), x\t),y^{t), 0) 

in D/G is a solution of the vakonomic equations for the constrained system (/, D/G) 
if and only if 

( d Bl f)l f)l 
d r)l 81 
81 

= - y'^eL)^ + (iVf„ - y^eDp^- 

Observe that if the basis of sections of D/G, {ci, Ca}, is orthonormal then 



and the vakonomic equations are: 

' = -T.(^'^% + y'l^Dy"" - {i'n, + yVrjpo, 

a 

c 
h 

Of course, it is also possible to study the abnormal solutions taking the lagrangian 
I = 0. 

Example 5.11. (See [19j[31]/ As a simple but illustrative example, consider 

with the distribution D = kerw where uj is the Martinet 1-forni: u — dx^ dx"^ . 

2 

Consider the vector fields generating D: 

8 {x^f 3 



3 



8x^ 2 3x-^ ' dx^ ' 
and the bundle metric { , )d which makes both vector fields orthonormal. 
Take now the action by translations: 

]R2 X — > ]R3 
{{a,b),{x^,x'^,x^)) I — > {x'^,x'^ + a,x^ + b) . 



We have a principal bundle structure tt 
M^-invariant. 



= K being both Z? and ( , )d 
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d d Ix^^'^ 3 
Observe that Vi: = span {^-j, ^~3^' ^ ^ -DHFtt = span {Vi = H ^ — J' 

d 

and = span {-;:r~r}- 
ax^ 

On the Atiyah algebroid TM^/M^ ^ M x R'"^ R we have the induced basis of 
sections : M ^ R x R^, z = 1, 2, 3: 

ei{x) = (x; (1,0,0)), 

62(2;) = (x;(0,l,y)) 

63(0:) = (a:; (0,0,1)) 

which induces coordinates {x,x,y^,y'^). Then, Z?/R^ is characterized by the equa- 
tion = 0. 

Observe that 

[61,62] = xea, [61,63] =0, [62,63] =0, 

and 

l{x;x,y') = laxf + {y'r). 
Consequently the vakonomic equations are: 

P2 = 

X = -y^xp2 
y^ = xxp2 . 

That is, P2 — k constant, and the vakonomic equations are precisely, 

X = —ky^x 
y^ = kxx . 

Since ^{x^ + (y^)^) is a constant of motion then we can take x{t) — rsm9{t) and 
y^{t) — rcosd{t) where 9{t) verifies the equation of pendulum 9{t) = —kr sm9{t) 
(see also Example l5.12p . 

5.3.5. Optimal Control on Lie algebroids as vakonomic systems. (See jl0[|29| ). Let 
T : E ^ Q he a, Lie algebroid and C be a manifold fibred over the state manifold 
TT : C Q. We also consider a section a : C ^ E along vr and an index function 

i-.c^m.. 

One important case happens when the section a : C ^ E along tt is an embed- 
ding, in this case, the image M = a{C) is a submanifold of E. Moreover, since 
a : C — M is a diffeomorphism, we can define L : M — > R by L = Z o . In 
conclusion, it is equivalent to analyze the optimal control defined by [l, a) (applying 
the Pontryaguin maximum principle) that to study the vakonomic problem on the 
Lie algebroid t : E Q defined by (L, M). 

More generally (without assuming the embedding condition), we can construct 
the prolongation : 7^C — > C of the Lie algebroid t : E ^ Q over the smooth 
map TT ; C ^ Q, that is 

T^C = {(6, Xp) e X TpC I pie) = T7:{Xp)}. 

Moreover, we have the constraint submanifold M characterized by 

M = {{e,Xp)e7^C\a{p)^e} 

and the lagrangian function L : 7^C — > R given hy L = Iot'^ . This is the vakonomic 
system associated with the optimal control system. If i? = TQ is the tangent bundle 
of the state space Q, it is not difficult to show that the prolongation of TQ along 
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TT : C —>■ Q is just the tangent bundle TC. Under this isomorphism, the constraint 
submanifold is 

M ^{X eTC\ Tt:{X) = cr{Tc{X))}. 
Thus, we recover the construction in |10j . Section 4. 

Example 5.12. Consider the following mechanical problem [51 HOI HI]- A (ho- 
mogeneous) sphere of radius r = 1, mass m and inertia about any axis fc^, rolls 
without sliding on a horizontal table which rotates with constant angular velocity 
about the x^-axis. The coordinates of the point of contact of the sphere with the 
plane are (x^,x'^). The configuration space of the sphere is Q = x 5*0(3) and 
the Lagrangian of the system corresponds to the kinetic energy 

K{x^,x'^;x^,x'^,ujxi,uJx^,^x^) = -^^f^ii^f + ^(i^)^ + mfc2(w^i + tj^2 + ^^la)), 

where {uxIjUJx^t'-^x'') are the components of the angular velocity of the sphere. 

Since the ball is rolling without sliding on a rotating table then the system is 
subjected to the afhne constraints: 

X^ — LOr^l = —Vlx^, 
x^ + Wj-i = Q,x^ , 

where VI is constant. Moreover, it is clear that Q ^M? x SO{3) is the total space 
of a trivial principal S'0(3)-bundle over and the bundle projection : Q ^ 
is just the canonical projection on the first factor. Therefore, we may consider the 
corresponding Atiyah algebroid TQ/SO{3) over M^. 

Since the Atiyah algebroid TQ/ SO{3) is isomorphic to the product manifold 
TM^ X so(3) ^ TR2 X R^, then a section of TQ/SO{3) ^ TR^ x R^ ^ is a 
pair {X,u), where X is a vector field on R^ and u : R^ ^ R^^ is a smooth map. 
Therefore, a global basis of sections of TR^ x R^ ^ R^ is 

d d 

ei = (^>0), 62 = (^,0), 63 - (0,£;i), 64 = (0,^2), 65 = (0,^3). 

where {Ei, E2, E^} is the canonical basis on R'^. 

The anchor map p : TR^ x R^ — > TR^ is the projection over the first factor 
and if |-, •] is the Lie bracket on the space Sec{TQ / SO{3)) then the only non-zero 
fundamental Lie brackets are 

[es, 64] = 65, |e4, 65] = 63, |e5, 63] = 64. 

It is clear that the Lagrangian and the nonholonomic constraints are defined 
on the Atiyah algebroid TQ / S0{3) (since the system is S'0(3)-invariant). In fact, 
we have a nonholonomic system on the Atiyah algebroid TQ / SO{2>) = TR^ x M^. 
This kind of systems was recently analyzed by J. Cortes et al [S] (in particular, this 
example was carefully studied). 

After some computations the equations of motion for this nonholonomic system 
are precisely 

X^ — LUx^ = —ftx^, 1 

x'^+ujx^ = ^x^, > (5-14) 
^3.3 = 6 J 

where c is a constant, together with 
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Now, we pass to an optimization problem. Assume full control over the motion 
of the center of the ball (the shape variables) and consider the cost function 

and the following optimal control problem: 

PLATE-BALL PROBLEM 2lj. Given points go, 91 e Q, find an optimal 
control curve {x^{t), x'^{t)) on the reduced space that steer the system from to qi, 
minimizes ^ {{x^)'^ + i^'^)'^) dt, subject to the constraints defined by Equations 

Observe that the Plate-Ball problem is equivalent to the optimal control problem 
given by the section : x ^2 ^ ^ ^^^^j^g j^^2 x R3 ^ R2 ^^^^^ 

(j{x^ , a;2; , u^) — (x^ , a;2; ^ . —u^ + flx^ , + Qx^ , c) . 

and index function l{x^ ,x'^;u^ jU^) — i((M^)2 + (^2)2). Since a is obviously an 
embedding this is equivalent to the proposed Plate-Ball problem. 

A necessary condition for optimality of the Plate-Ball problem is given by the 
corresponding vakonomic equations. For coherence with the notation introduced 
along this paper, denote by 

= y2 ^ ^2^ ^3 ^ ^ ^4 ^ ^ ^5 ^ 

Therefore, the vakonomic problem is given by the Lagrangian L = ((2/^)^ + (j/^)^) 
and the submanifold M defined by the constraints: 

y3 ^ ^^x\x\y\y^) = -y^ + nx\ 
y4 ^ ■^\x\x^,y\y^) ^y^ +nx'^, 
= (x^ ^ x'^ , y^ , y'^) = c 

After simple computations we obtain that the vakonomic equations are: 

P3 = cp4 - {y^ + nx"^) P5 

P4 = ~cp3 - (y2 - nx^^ P5 

P5 = (y^ + ^x^) Pi + {y^ - ^x^) Pi 

J^ijj'-Pi) - -^Pi 



y^ = x^ , y^ = ±2 



The system is obviously regular since the matrix 

f d^L A 92^a / 1 



1 dyaQyb « dyaQyb / V ^ ^ 

\ a=3 y y / ,,<2 \ 

is non-singular. Therefore, there exists a unique solution of the vakonomic equations 
on Wi, determined by the conditions: 

dL d^" _ 1 
dy^ dy^ 
dL 8^°' _ 2 
dy'^ 9?/2 

Therefore, it follows that {x^^ ,x'^ ,pi,p2,Pi,Pi,P5) are local coordinates on Wi (or, 
from Corollary 5.9, on E* , if you prefer). 



= 7rT~P<^jrT = y -Pi 



P2 = -^-Pa^TT^y +P3 
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Moreover, on Wi we have a well-defined Poisson bracket { , }(L,Af) whose non- 
vanishing terms are: 

{x^,Pl}(L,M) = 1 , {x'^,P2}(L,M) = 1 
{P3,P4}{L,M) = -P5 , {P3,P5}{L,M) = P4 , {P4,P5}(L,A/) = -PS 

In these coordinates the Hamiltonian Hwi is; 

i?i4'i(a;\a;^,Pi,P2,P3,P4,P5) = ^{pi + P4f + ^{P2 -P3f + cp^+P3^x^ -{^PiVlx^ 
and the vakonomic equations are 



Pi 


= {Pli Hwi}(L,M) 


= -^P3 


P2 


= {P2,Hwi}(L,M) 


— — rip4 


P3 


= {P3i Hwi}(L,M) 


= cp4 - (pi +P4 + ilx'^) P5 


Pi 


= {Pii Hwi}(L,M) 


= -CP3 - {P2 -P3- ^x'^) Pb 


P5 


= {P5,Hwi}(LM) 


= {pi +PA + ^X^) P3 + {P2 - 




= {x^ , Hwi}{L,M) 


= Pl+PA 




= {x'^,Hwi}{L,M) 


= P2-P3 



One of the most studied cases is when c = (the sphere is rolled that its angular 
velocity is always parallel to the horizontal plane) and = (not rotation of the 
plane). In this case the equations reduced to 

P3 = -y^Pb ,P4 = -y'^Pb , Pb = y^P3 + y'^PA 

1 -1 2 -2 

y = X ,y = X . 
From these equations it is easy to deduce that 

j^iiy'r + iy'f) = o 

Therefore y'^{t) = cos0{t) and y'^{t) = sme{t), when 1 = ^/iy^^^+Jy^ . More- 
over, it is easy to deduce that p^ = 9. 

Now, since y^ = p4 + fci and y^ = —ps + k2 with fci,A:2 G R constants, taking 
ki = r cos (fi and k2 — r sin ip then from Equation p^ = y^P3 + y^Pi we deduce that 

9 — r cos 9 sin ip — r sin 9 cos (p — — r sin(0 — (p) 

that is, the angle 9 satisfies the equation of pendulum, while the coordinates of the 
contact point satisfy the ODEs: x^ = cos9 and x^ = sm9. The remarkable result 
is that these equations say that the contact point of the sphere rolling optimally 
traces an Euler elastica (see [H]). 

6. Conclusions and future work 

We have developed a general geometrical setting for constrained mechanical sys- 
tems in the context of Lie algebroids. We list the main results obtained in this 
paper: 

• We develop a constraint algorithm for presymplectic Lie algebroids and 
discuss the reduction of presymplectic Lie algebroids. 
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• For a singular Lagrangian function on a Lie algebroid, we look for solu- 
tions of the corresponding dynamical equation, applying the previously 
introduced constraint algorithm. In addition, we find the submanifolds 
where the solution is a SODE (a second order differential equation). The 
theory is illustrated with an example. 

• We study vakonomic mechanics on Lie algebroids. As for singular La- 
grangian systems, we deduce the vakonomic equations by means of a cons- 
traint algorithm. We define the vakonomic bracket in this setting. Fur- 
thermore, the variational point of view and some explicit examples are 
discussed. We remark that we just look for the so-called normal solutions. 
We postpone for a future work a detailed variational analysis of vakonomic 
mechanics on Lie algebroids, including abnormal solutions... We will illus- 
trate our theory with an example related with Optimal Control Theory, 
but our approach may cover other interesting examples (see [21 [15]) and 
admits, in a standard way, variational discretizations (see [21 117|). 
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